Reducing decoherence of the confined exciton state in a quantum dot by 

pulse-sequence control 
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We study the phonon-induced dephasing of the exciton state in a quantum dot excited by a 
sequence of uhra-short pulses. We show that the multiple-pulse control leads to a considerable 
improvement of the coherence of the optically excited state. For a fixed control time window, the 
optimized pulsed control often leads to a higher degree of coherence than the control by a smooth 
single Gaussian pulse. The reduction of dephasing is considerable already for 2-3 pulses. 



I. INTRODUCTION 

Modern semiconductor technology offers a range of 
ways to produce artificial systems where a small num- 
ber of charges, confined in all three spatial dimensions, 
form atomic-like systems. These quantum dotsA^^ (QDs) 
allow one to optically control the quantum states in a 
way typical for atomic system a'^i^i^i^i'^'i^i^i-^" , but they are 
still semiconductor charge devices, opening the possibil- 
ity of integration into future nanoelectronic and quantum 
electronic devicesii^i^. 

The recent progress in ultrafast spectroscopy of semi- 
conductor systems^^ made it possible to control and 
probe the quantum states of carriers confined in a QD 
on femtosecond time scales. However, any interaction 
of the confined carriers with the external driving fields 
not only leads to the desired quantum transitions but 
also can induce some unwanted ones. If both the de- 
sired and unwanted transitions have a discrete nature 
(e.g., the exciton vs. biexciton transition in a QD) the 
latter may be suppressed by a suitable choice of control 
pulsesi^ii^. However, since the QDs are embedded in the 
macroscopic crystal, an important role is played here by 
phonon-assisted transitions which inherit the continuous 
nature from the phonon states. The presence of a contin- 
uum in the system state obtained after optical excitation 
leads to irreversibility of the subsequent kinetics and to 
dephasing of the quantum statesiSiiLi^iiSiSS even in the 
absence of real transitions (energy relaxation). In the 
latter case these processes are known as pure dephasing. 

Effects of pure dephasing may be observed in mea- 
surements of the optical polarizationiiiSiiS or exciton 
occupation in a QDi2i2ii2^. The qualitative agreement 
between theoretical predictions based on the phonon- 
induced pure dephasing modeliS and four-wave-mixing 
data22, confirms that these phonon-related processes 
dominate the dephasing in the sub-picosecond resonant 
driving regime. Recently, it has been shown that the 
initial decoherence in QDs with strong electronic con- 
finement as described by the pure dephasing model 
is also in quantitative agreement with corresponding 
measurements^^. It turns out that even at low temper- 
atures (~ 5 K) typically about 20 % of the coherence is 
lost already in the initial phase of the dynamics within 
a few picoseconds. Obviously, this rapid decoherence is 



not negligible and can potentially become an obstacle for 
many proposed device applicationsS^*^. 

The question arises whether it is possible to reduce the 
phonon-assisted component of the final state by some 
pulse shaping techniques, e.g., by using sequences of 
pulses, in analogy to the discrete transition casei^ti^. It is 
clear that the continuum nature of these unwanted tran- 
sitions makes this task much more demanding. One obvi- 
ous way of excluding the phonon-assisted transitions is to 
make the control pulse spectrally selective by extending 
the pulse durations. This is, however, disadvantageous 
when other decoherence processes are present that de- 
crease the coherence properties of the final state if the 
control operation takes too long and lead to a trade-off 
situation^^. For practical purposes it is therefore essen- 
tial to minimize the phonon-related pure dephasing with 
the constraint that the time window for the control action 
is limited. 

In this paper we propose a solution to this optimiza- 
tion problem in terms of a control sequence composed of 
a series of ultrashort (but finite) pulses. The advantage 
of this approach is that it is simple both in theoreti- 
cal formulation and experimental realization, as it does 
not require the technical ability to obtain arbitrary pulse 
shapes. Moreover, sequences of temporally broadened 
pulses have simple spectral features, even in the nonlin- 
ear regime, which lead to a transparent interpretation of 
the resulting dephasing effect. 

On the other hand, in spite of the apparently limited 
class of pulse shapes considered, the proposed approach 
turns out to be unexpectedly efficient. Already for a few 
pulses the degree of decoherence falls down considerably, 
while a further increase of the pulse number leads only to 
negligible improvement. We shall show also that achiev- 
ing by a single Gaussian pulse the same dephasing as that 
resulting from a sequence of a few pulses often requires 
much longer control times. Thus, the quality of coher- 
ent optical control of the excitonic system can indeed be 
increased by simple means, using series of phase-locked 
laser pulses, without the need to generate pulses of arbi- 
trary shape. 

The paper is organized as follows. In the next section 
we present the model of the carrier-phonon system. In 
Sec. mil we discuss the limiting case of infinitely short 
pulses which can be treated analytically. In order to 
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demonstrate the reduction of decoherence by multiple- 
pulse control we compare the coherence losses occurring 
after a single pulse excitation with a given pulse area with 
those occurring after a pulse sequence with the same total 
pulse area. Next, in Sec. IIVI we present a perturbative 
method for determining the coherence loss for arbitrarily 
shaped pulses and apply it to the case of multiple-pulse 
sequences. The loss of coherence after a series of driv- 
ing pulses with fixed finite pulse durations and optimized 
intensities is discussed in Sec. |V] Physically, the dura- 
tion of the pulses is limited to the lower end by the re- 
quirement to spectrally avoid higher energetic excitations 
which would be a source of further undesired dynamics 
and which are not accounted for in our model. In Sec. EH 
we extend our material model and subsequently include 
LO phonons and higher excited electronic states. Explic- 
itly including these excitations in our optimization pro- 
cedure allows us to use shorter pulses as building blocks 
of our pulse sequences. The undesired dynamics result- 
ing from excitations of these higher energetic states is 
now avoided by optimizing the pulse sequence. Finally 
Sec. IVIII concludes the paper. 



Effective dielectric constant 
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62.6 


Longitudinal sound speed 


Cl 


5600 m/s 


LO phonon frequency 


n 


54 ps-i 


Deformation potential 






electrons 




-8.0 eV 


holes 




LO eV 


Crystal density 


p 


5360 kg/m^ 


Wave function widths: 






electron in-plane 


le 


4.4 nm 


hole in-plane 




3.6 nm 


electron/hole in 2-direction 




LO nm 



TABLE L The GaAs material parameters and QD system 
parameters used in the calculations (partly after Refs. 



phonons will be considered, with the corresponding cou- 
pling constants given by^— 



Fk = 



I hk 

2pvc\ 



[<7,T,{k) - aMk)] (LA), (2) 



and 



II. THE MODEL 

The environment-induced process of loosing the phase 
relations between different components of quantum su- 
perpositions may be observed experimentally in a num- 
ber of ways. One of its manifestations is the decay of 
the optical polarization excited by a certain laser pulse 
or pulse sequence resonantly coupled to the ground-state 
excitonic transition in a quantum dot. The impact of this 
dephasing process was studied in the case of lineariiSiii 
and nonlineariiSiiS optical experiments. 

Here, we consider the usual model of carrier-phonon 
interaction involving the two states of the carrier sub- 
system: the empty dot, |0), and the one exciton state, 
1), the latter being coupled to phonon modes. The wave 
functions will be modeled by anisotropic Gaussians. We 
assume a simple product (Hartree) form of the exciton 
wave function with the hole component shrunk in-plane 
due to Coulomb interactions (this is consistent with pre- 
vious results that were obtained by full numerical diago- 
nalization methodsSS). 

Thus, in the rotating wave approximation, the Hamil- 
tonian reads (in the rotating frame) 



H = -f{tm{i\ 



|l)(0|) + |l)(l|^F;,(&,+6Lj, (1) 



where f(t) is the pulse envelope (which will be assumed 
real), bk are the phonon annihilation operators for the 
mode k (including implicitly also the branch index) and 
Fk are phonon coupling constants for the ground exci- 
ton state. In the present paper the deformation potential 
coupling to the longitudinal acoustical (LA) phonons and 
the Frohlich coupling to the longitudinal optical (LO) 



-^[T,{k)-Tuik)] (LO), (3) 



where p is the crystal density, ci is the longitudinal 
sound speed, v is the normalization volume of the phonon 
modes, e is the elementary charge, eo is the vacuum di- 
electric constant, e = (l/coo — l/cs)"^ is the effective 
dielectric constant, where eg and Eoo are the static and 
high-frequency dielectric constants, respectively, Q is the 
LO phonon frequency and 



J^o.hik) 



d^rvl/*jr)e*'-vl',,h(r). 



where ^'e,h('") are the electron and hole wave functions. 
Since the electron and hole are assumed to overlap, the 
piezoelectric coupling will be neglected^i. The parame- 
ters used in the computations are shown in Tab. ^ 



III. POLARIZATION DEPHASING AFTER A 
SEQUENCE OF ULTRASHORT PULSES 

The natural starting point for the discussion of the 
multiple-pulse control of the optical polarization in a QD 
seems to be the limit of ultrashort pulses. In this case, 
each pulse in the sequence acts on a time scale much 
shorter than the lattice response times. The subsequent 
evolution corresponds to a relaxation of the lattice to a 
new equilibrium that is determined by the optically cre- 
ated confined charge distribution. During this process 
the exciton occupation cannot change but the carrier- 
phonon correlations that appear in the system partly de- 
stroy the exciton coherence. The resulting decay of the 
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FIG. 1: The evolution of the optical polarization excited with 
one or four infinitely short pulses. Solid: exact result, dashed: 
perturbative result. 



optical polarization may be described analytically, both 
for a single pulse^ and for a sequence of pulsesiS,. 

It has been shown experimentally-^'-ii that the optical 
polarization excited by an ultra-short pulse undergoes an 
initial dephasing on a picosecond time scale, followed by 
a much slower exponential decay. Model calculations re- 
veal that the initial decay in these systems can be mostly 
attributed to LA phonons^^'^'*. In general, decoherence 
refers to the loss of definite phase relations between vari- 
ous components of a quantum superposition. In our case 
these components are the crystal ground state and the 
QD exciton state. Most easily affected by decoherence is 
the superposition state with equal amplitudes, which is 
reached from the ground state by an a = 7r/2 rotation on 
the Bloch sphere. In the following we will concentrate on 
minimizing the decoherence occurring after such a 7r/2 
rotation. Other optimization goals have been considered 
and treated with other strategiea^^. 

Analytical results for LA phonon-induced pure dephas- 
ing obtained along the lines described in Ref . are plot- 
ted in Fig. Shown is the evolution after a single 7r/2 
pulse and after a series of four identical, in-phase, equidis- 
tant 7r/8 pulses. The total length of the series is ttot = 1 
ps and the last pulse arrives at t = 0. The spectral 
cut-off resulting in real experiments from the finite pulse 
duration is introduced in this (5-pulse model by includ- 
ing only the acoustic phonon response, which amounts 
to the assumption that the actual pulse duration is still 
longer than LO phonon periods and also long enough to 
spectrally avoid higher excitations. 

Compared to the single-pulse excitation, the initial po- 
larization in the 4-pulse case is reduced due to dephasing 
in between the pulses. However, the final level of po- 
larization is higher. This means that the multiple-pulse 
excitation leads to a lower contribution of the phonon- 
assisted transitions in the final state. At zero tempera- 
ture the loss of coherence compared to the single-pulse 
initial value is reduced from about 12% to about 9%, i.e., 
we obtain an improvement of 25%. In the remaining part 
of the paper we will explore the possibilities of reducing 
the degree of dephasing, going beyond the (5-pulse limit 
and, finally, extending the modeling to higher-energetic 
excitations. 



IV. PHONON-INDUCED PERTURBATION OF 
THE COHERENT DYNAMICS 

It may be expected that the reduction of dephasing 
described in the previous Section may be extended by 
increasing the number of pulses and optimizing their am- 
plitudes. However, the analytical formulas become highly 
involved for growing numbers of pulses. In any case the 
spectral width of a sequence of 6 pulses will be infinite 
(although, possibly, the spectra may get non-trivially 
shaped). Here, it should be recalled that the spectral 
selectivity of real experiments has been accounted for in 
this approach by concentrating on the lattice response in 
the low-frequency regime and disregarding in the mate- 
rial model all higher energetic excitations. A more ac- 
curate treatment of the spectral selectivity is, of course, 
to introduce a spectral cut-off. This requires using finite- 
length pulses; a case for which within our model no exact 
solution of the dynamics is known. 

In order to deal with these issues and to continue 
and extend our discussion without restricting it to lin- 
ear optical effects we resort to a scheme which is non- 
perturbative in the electromagnetic field but includes the 
phonon effects only in the leading order. 

To calculate the phonon-related perturbation we ex- 
pand the formal expression for the density matrix of 
the total system (in the interaction picture) at the fi- 
nal time t up to the second order in the phonon coupling 
constants^. 



(4) 



^n■ Jti 

- ^ £ dT J^^ dr' [iJint (r) , [i/int {t'),q{U)]1 

where Q{t) and Hi-^tir) are, respectively, the density ma- 
trix of the total system and the carrier-phonon interac- 
tion Hamiltonian in the interaction picture with respect 
to the laser-induced evolution (which, in the case of res- 
onant coupling, is known exactly) and is the initial 
time. We assume an uncorrelated initial state with the 
lattice in thermal equilibrium and the exciton subsystem 
in the ground state. This procedure does not involve any 
Markovian approximations; the truncation at the second 
order is valid as long as the correction to the density 
matrix remains small. 

After tracing out the lattice degrees of freedom and 
setting a = 7r/2 the polarization at the time t, which is 
proportional to goi (i) i may be written as 



\P{t)? = \Po\ 



1 _ / du:^ \K{u:) 



where 



K{uj) = / dre'"^^— sin$(T), 
dr 



$(t) 



dr'fir'). 



, (5) 



(6) 
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The upper limit of integration in Eq. JBJ has been ex- 
tended to infinity under the assumption that at time t 
the pulse has already been completely switched off (the 
details of the formalism are described in Ref. l34l). Since 



the function R{uj)/uP' is regular dX uj — Q, the asymptotic 
limit of the polarization at t oo is (cf. Rcf. 17) 



|PMI 



l^^ol 



1 



1 



(7) 



The phonon spectral densities R{uj) characterize the 
spectral properties of the lattice and are given by 



k 



(8) 



- ujk) + 5{uj + ujk)] 



where n-aiuj) is the Bose distribution function. 

It should be noted that part of the polarization drop 
can be attributed to the deformation polaron formation 
around the exciton (phonon dressing) which is not nec- 
essarily irreversible. In fact, a perfectly adiabatic, thus 
reversible, passage to the exciton-no-exciton superposi- 
tion corresponds to vanishing K{oj) and, although no ir- 
reversible dynamics takes place, there is still a relative 
polarization reduction 



R{lo) 



which corresponds to 1/2 of the asymptotic value after a 
5 pulse. 

Since the spectral function K{u}) [Eq. ©] is a nonlin- 
ear functional of the pulse envelope f{t), the only way 
of calculating it is, in general, to numerically perform 
the Fourier integral. However, some insight in its struc- 
ture is possible in the case of a series of non-overlapping 
control pulses. Let us denote the nth pulse envelope 
(n — 1, . . . ,7V) by /„(< — i„), where t„ is the time at 
which the pulse is applied, and the value of $(i) after 
the nth pulse by a„, with ao = and = a. Then, 
during the nth pulse one has 



$(r) =a„_i + <l„(T-t„), 



where 



fnir')dT' 



Substituting this into Eq. ^ one finds 



K{co) = [F{uj)~F*i-u)], 

N 



with 



F{u;) = ^F„(^)e^("-i+'^ 



(9) 
(10) 
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FIG. 2: Spectral characteristics of the phonon reservoir and 
of the driven dynamics. Solid line: The spectral density 
R{uj)/llP at T = K for LA phonons. Dashed line: the 
nonlinear spectrum |_ftr(tj)p for a sequence of identical short 
pulses for a — tx/2 and pulse numbers and total durations 
as shown. Dash-dotted line: the envelope calculated as the 
linear spectrum of a single pulse (up to an appropriate fac- 
tor). In the upper figures also the nonlinear spectrum for a 
sequence of (5- pulses is shown (dotted). 



Since, by definition, the function <I'n(i) changes only 
around t = 0, the transform Fn{Lij) is a smooth enve- 
lope of width inversely proportional to the duration of 
the pulse fn{t)- If the rotation angle per pulse is small 
this envelope function is close to the linear Fourier trans- 
form of /n(t) (up to a factor). 

In the limit of infinitely short pulses the transform 
Fn{Lu) becomes infinitely broad and one has 



F„(w) = F„(0) = e 



1. 



(11) 



If, in addition, the pulses are equally spaced, i„ = nAt, 
then the resulting spectral function F{lu) is periodic in 
the frequency domain with the period 27r/Ai. This limit 
corresponds to the exact solution obtained previously and 
allows one to compare both results and estimate the error 
introduced by the perturbative approach. As can be seen 
in Fig. ^ the perturbative results not only reproduce 
qualitatively all the essential features of the exact result 
but are also reasonably close to the exact results, at least 
at low temperatures. 

Temporally broadening the pulses results in the ap- 
pearance of the envelopes Fn{uj) [Eq. (|10|l ] which break 
the periodicity and eliminate the high-frequency part of 
the nonlinear spectrum. For a series of identical, weakly 
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overlapping, equally spaced pulses one has Fn{uj) = 
F{uj), an = na/N, and = {n — l)f/N. In order to 
simplify the formulas, we define here T = Nttot/{N — 1), 
where ttot is the total length of the pulse sequence mea- 
sured from the center of the first pulse to the center of the 
last. The summation in Eq. may be easily performed 
analytically and the result for a = 7r/2 is 



(12) 



sin' 



sin2^(^ + ^) 



sin^i(--i.) 



sm 



2T' 



This is formally identical to the sum of two diffraction 
patterns for equally spaced optical grids of length T, 
modulated by the envelope F(uj) which suppresses higher 
order "diffraction maxima", and shifted oppositely by 
±7r/(2T). This behavior is illustrated in Fig. 12 where 
we have chosen pulses with a duration of Tp — 100 fs 
(full width at half maximum, FWHM). As seen in the 
figure, this choice assures that the relevant spectral func- 
tion |A'(tj)p has a spectral cutoff below the LO phonon 
frequency = 54 ps""'^, consistent with the neglect of 
this excitation in our model. 

In this case of identical, weakly overlapping, equally 
spaced pulses, for increasing pulse number the period of 
the spectrum in the frequency domain increases, i.e. the 
higher maxima are shifted towards higher frequency and 
decrease due to the envelope. However, the width of the 
central maximum is determined by the total duration of 
the sequence. From Eq. H12() it can be seen that when 
taking T as a measure of the sequence duration then 
this width is essentially insensitive to the pulse number. 
Similarly, the value at w = is fixed by the gate angle 
and is always |i^(0)|^ = sina = 1. Modifying the total 
duration of the sequence leads to proportional scaling of 
the w-dependence of the nonhnear spectrum |i4r(w)l^ (the 
"diffraction pattern") under the envelope. The latter is 
fully determined by the single pulse spectrum which, for 
a < 7r/2 is very close to the linear one even for two pulses 
(see Fig.ig. 

The independent control over the nonlinear spectrum 
for high frequencies (by pulse width) and in the LA fre- 
quency domain (by the sequence length and pulse num- 
ber) opens the possibility of optimizing the coherence by 
reducing the overlap in Eq. , alternative to using long 
smooth pulses. In the next Section we study this issue, 
releasing the restriction to equal pulses in order to find 
optimized pulse sequences. 



V. REDUCTION OF DEPHASING BY 
OPTIMIZED PULSE SEQUENCES 

The aim of this Section is to quantify the reduction 
of the LA-phonon induced dephasing of the optical po- 
larization by using optimized sequences of pulses with a 
fixed finite pulse length. As discussed in the previous 
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FIG. 3: The relative polarization drop (with respect to the 
unperturbed case) for the a = 7r/2 excitation performed by 
sequences of A'^ = 2, 3, 4 equidistant narrow pulses with Tp = 
100 fs within the time interval ttot at T = K for optimized 
pulse intensities (solid) and for identical pulses (dashed). 



Section, by using pulses of finite duration we impose an 
envelope on the nonlinear spectra of the controlled dy- 
namics and thereby assure that higher-frequency spectral 
features, including LO phonons, higher exciton states and 
confined two-pair excitations, have negligible impact. 

In real systems, apart from the pure dephasing result- 
ing from the phonon response to the evolution of the 
exciton state there is another contribution to the loss of 
fidelity related to exponential decay processes (finite ex- 
citon lifetime in our specific case)^'^''^^. The effect of the 
latter on the polarization magnitude grows with the du- 
ration of the pulse sequence. Therefore, a question of 
practical importance is to what extent the system coher- 
ence may be preserved against the pure dephasing during 
an operation performed by a sequence of a fixed number 
of pulses within a definite time interval that is limited by 
the exponential decoherence. 

By numerically calculating the nonlinear spectra K{ll!) 
[Eq. ©J for sequences of equidistant narrow (but not 
necessarily non-overlapping) pulses and maximizing the 
asymptotic polarization magnitude [Eq. (jT))] with respect 
to the pulse amplitudes for a fixed total rotation angle 
a — TT/2we found optimized^ pulse sequences for a given 
total duration of the sequence itot and the resulting min- 
imal polarization drop 



d = 



W-\Pi^)\' 



The result is shown in Fig. |21 

As expected, it turns out that the degree of dephasing 
typically tends to decrease as the time interval becomes 
longer. Optimizing the pulse sequence leads to some fur- 
ther reduction of the polarization drop compared to the 
case of equal pulses, although this additional gain is rela- 
tively small. An interesting feature is the slight increase 
of the decoherence drop for times ttot ?J 1 ps in the N = 2 
case. This may be understood with the help of Fig. [21 
by noting that for 27r/ttot — ^^LA, where wla is the LA 
spectral density cut-off, the maximum of R{uj)/uj'^ coin- 
cides with the minimum of [if (w)!^ while for longer times 
the latter is quickly oscillating and averages to its mean 
value. 
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FIG. 4: The optimized sequences of = 2 (upper) and A'' = 4 
(lower) narrow pulses (rp — 100 fs) compared to the Gaussian 
pulses (dashed, arbitrary scale) yielding the same asymptotic 
polarization (at T = K). 



For itot ^ Tp the duration of the pulse sequence is ac- 
tually determined by Tp and ftot becomes the shortest 
relevant time scale. Thus, the high frequency part of the 
pulse spectrum is governed by itot /N in this case and the 
simple discussion presented above for non-overlapping 
pulses is not valid. In this case, the nonlinear spectra 
extend to much larger frequencies and the overall error 
drops down due to strong modulations of the resulting 
pulse shape. Although the spectral cut-off cannot be as- 
sured a priori in this case, it is still possible to control the 
consistency of the calculation by inspection of the non- 
linear spectrum corresponding to the resulting optimized 
pulse sequence. In fact, it turns out that for very short 
ttot the error values shown in Fig.|31do not correspond to 
pulse sequences with properly restricted spectra. More- 
over, due to the very complex structure of the error as a 
function of the pulse intensities it is not guaranteed that 
the optimized value found by our procedure in the strong 
overlap regime corresponds to the global minimum. For 
weak and moderate overlap the optimized sequences are 
likely to correspond to the global minimum. 

It is interesting to compare the coherence loss for an 
excitation performed by a sequence of pulses and by a 
single Gaussian pulse. In Fig. ^ we compare the opti- 
mized pulse sequence with a Gaussian pulse performing 
the same a = n/2 rotation with the same asymptotic po- 
larization drop as obtained from the sequence. It turns 
out that already for = 2 the pulse sequence is shorter 
than the corresponding Gaussian pulse as long as itot ^ 2 
ps. For very short pulse durations the Gaussian may 
even be twice as long as the pulse sequence. For itot ^ 2 
ps, the smooth Gaussian control can no longer be sig- 
nificantly outperformed by a pulse sequence. However, 
as the number of pulses is increased the pulsed control 
becomes more efficient over a wider range of sequence 
durations. 

Another interesting aspect worthwhile to note is that 
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FIG. 5: The polarization drop as a function of the number 
of pulses N for the total durations ttot = 1 ps (squares) and 
= 5 ps (circles) at T = K (left) and T = 10 K (right). 
Full symbols: optimized pulse intensities, empty symbols: 
equal pulses. Here Tp = 100 fs. 



in the infinitely short pulse limit the optimal sequence 
of A^ = 2 pulses does not depend on the structure of the 
spectral density of the reservoir and for a — 7r/2 it always 
corresponds to the ratio of pulse areas of 1:2 (see the 
Appendix). As can be seen in Fig.^ this universal area 
ratio of the optimal 2-pulse sequence holds approximately 
also for pulses of finite duration. 

By increasing the number of pulses, the coherence loss 
may be further reduced (Fig. [S]) . We find that at both 
temperatures the polarization drop can be reduced by 
about 40%. Note that even in the adiabatic limit, 50% 
of the polarization drop is present. This implies that 
the irreversible effect is reduced by about 80%. The re- 
duction strongly depends on the total duration but in all 
cases (for ttot in the picosecond range) most of it is gained 
already for a few pulses. For short itot and for a control 
sequence consisting of equal pulses there is an optimal 
number of pulses beyond which the error grows again. 
This results directly from Eq. H12|) : for higher A^ the 
nonlinear spectrum decreases more slowly around uj = 
which determines the resulting overlap with R{uj) when 
the width of the "zeroth order" maximum is comparable 
to the LA cut-off frequency. This effect disappears if one 
allows variable pulse intensities. From Fig. |S1 it is evi- 
dent that for sufficiently long total durations ttot of the 
pulse sequence the polarization drop S comes close to its 
value in the adiabatic limit already when the sequence 
comprises only a few pulses. In this case, the optimaza- 
tion of the intensities yields only a marginal improvement 
as is seen from the curves corresponding to ttot = 5 ps 
in Fig. |31 In contrast, for pulse sequences of short du- 
ration ttot, 6 is well above its adiabatic value, but the 
improvement resulting from the intensity optimization is 
now substantial as illustrated by the ttot = f ps curves 
in Fig.El 

In order to see most clearly how the optimization of 
overlapping pulses leads to a decrease of decoherence 
let us compare the spectral features corresponding to 
equal and optimized pulses in the case of strong over- 
lap (Fig. 13. In the optimized strong modulation 
of the resulting pulse envelope is possible which leads 
to considerable reduction of the nonlinear spectrum in 
the LA sector. Such a modulation produces, however. 
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FIG. 6: Left: the optimized pulse envelope formed by the se- 
quence of Af = 4 overlapping narrow pulses (dashed: pulses 
with optimized amplitudes, dotted: equal pulses) for the se- 
quence duration ttot = 0.3 ps at T = K. Right: the cor- 
responding nonlinear spectra [/("(lj)!^ along with the phonon 
spectral density (solid). 

large spectral features beyond the LA cut-off, as can be 
clearly seen in Fig. El In order to provide the desirable 
spectral cutoff the pulse length must now be increased, 
compared to the case of non-overlapping pulses, or fur- 
ther constraints have to be introduced in the optimiza- 
tion procedure. 



VI. EXTENDING THE MODEL 

As discussed above, the pulse lengths used to control 
the quantum state must be long enough to assure a fre- 
quency cut-off that avoids in any case the undesired dy- 
namics that otherwise may result from high energetic ex- 
citations. This restriction may be inconvenient for prac- 
tical applications as it implies the use of long pulses as 
building blocks of the control sequence. Below we show 
how the optimization procedure may be extended such 
that the undesired dynamics is not avoided by a fre- 
quency cut-off but rather by suitably shaping the pulse 
sequence. To this end wc first need to extend our model 
and include explicitly the pertinent higher energetic ex- 
citations. 

The first obvious step is to include LO phonons which, 
for confined states, contribute a single Hne (due to their 
weak dispersion around fc = 0). In the plot of the po- 
larization drop against the duration (Fig. [7| left) small 
additional features appear periodically, both for equal as 
well as for optimized pulses (visible in the inset). They 
result from the fact that in the case of weakly overlap- 
ping pulses the positions of the maxima in the nonlinear 
spectrum [^(w)!^ are fixed by the total duration ttot and 
the total number of pulses, and the decoherence effect 
increases each time one of them overlaps with the LO 
phonon line (Fig. [7| right). In the system under dis- 
cussion this effect is rather weak because of the strong 
reduction of the LO phonon line due to charge cancella- 
tion. Note that extending the reservoir model allows us 
to reduce the pulse duration to 60 fs which now provides 
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FIG. 7: Left: the asymptotic polarization decay induced by 
LA and LO phonons for N = A pulses and rp = 60 fs; solid: 
optimized sequence, dashed: equal pulses. The inset shows 
part of the curves in magnification. Right: The nonlinear 
spectrum |A'(a;)p corresponding to two values of the pulse 
sequence duration ttot = 1.74 ps (strong overlap with the LO 
line; dashed) and ftot = 1.92 ps (weak overlap; dotted). Also 
shown are the spectral densities for LA (solid) and LO (bar, 
arbitrary scale) phonons. 

the spectral cutoff before the onset of higher-frequency 
reservoir features. 

The next higher excitations of the phonon bath are 
related to the possibility of phonon-assisted transitions 
to dark states^ (which may be viewed as an effect of 
nonadiabaticity of exciton dynamics with respect to LO 
lattice oscillations'^^). We extend the Hamiltonian ^ by 
adding terms describing the excited exciton states and 
their LO-phonon-assisted coupling to the ground state 
(other couplings do not contribute to the leading order), 

n>l n>l fc(LO) 

where £"„ are the energies of the excited states and the 
coupling constants i^^"' are analogous to those given by 
Eqs. Q and and may be found e.g. in Ref. 1201 
This coupling induces additional decoherence with the 
corresponding asymptotic polarization drop given by 

A|P(™)|= = 

-iP.i^|;/.„i^^(^(i.Hi^4), 

where 

f°° ■ d 1 
L{lu)= J dre-^_sin-$(r), (13) 

and i?^"''(a;) are defined as in Eq. JSJ but with the cou- 
pling constants F^"^ (see also Ref. 0). The optimal 
control, as measured by the minimal polarization drop, 
corresponds now to minimizing the sum of the overlap 
integral appearing in Eq. Q) and that given above. 

For the results shown in Fig. |S1 we include three opti- 
cally inactive excited states at 17, 32 and 33 meV above 
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FIG. 8: The pulse sequence (left) and the nonlinear spec- 
tra (right) corresponding to the optimal control in the model 
extended by adding LO-phonon-assisted excitation of a few 
lowest dark states (X1-X3). In the right panel, the solid and 
dashed lines show the spectral functions \K{uj)\^ and \L{u)\^ , 
respectively, while the vertical bars represent the spectral fea- 
tures related to phonon-assisted transitions to optically inac- 
tive states. 



the ground state. The optimization was done for a rel- 
atively high number of pulses, so that they overlap giv- 
ing rise to a single, strongly shaped envelope. The sim- 
ple "diffraction pattern" argument does not hold in this 
case and strongly nonlinear effects become important. As 
a result, additional flexibility appears in the shaping of 
the nonlinear spectra and the pulse optimization leads to 
spectral functions that to a large extent avoid the overlap 
with these discrete transitions as can be seen in Fig. |S| 



the resulting dephasing of the exciton states (measurable 
via the decay of the optical polarization). This is due to 
the fact that the whole continuum of phonon frequencies 
cannot be avoided simultaneously. Nevertheless, with an 
optimized choice of pulse intensities the degree of dephas- 
ing may be considerably reduced, both at zero and finite 
temperatures, to a degree depending on the allowed du- 
ration of the pulse sequence (as expected, longer control 
time windows admit larger reduction of dephasing). 

From the practical point of view, the results obtained 
in this paper are remarkable with two respects. First, 
in the picosecond range of control times (relevant for 
phonon-induced pure dephasing effects) using sequences 
of a few (2-4) pulses leads to higher final degree of coher- 
ence of the resulting state than that obtained by a single 
Gaussian pulse of the same total duration. This may 
be important since obtaining short series of phase-locked 
pulses is certainly much more feasible than generating an 
arbitrarily shaped pulse resulting from a full optimization 
procedure. 

Second, it is found that increasing the pulse number 
beyond 3 or 4 brings only small improvements, compared 
to that gained by replacing the single pulse by a series of 
2 or 3. This shows that the limited possibility of reducing 
the dephasing may be almost fully exploited already with 
limited control resources. The universality of the optimal 
sequence of two pulses further increases the usefulness of 
these results. 



VII. CONCLUSIONS 

We have shown that the coherence of the optically gen- 
erated superposition of exciton occupation states con- 
fined in a quantum dot may be considerably increased 
by using a sequence of narrow (but finite) pulses to drive 
the carrier system. By doing so, one modifies the nonlin- 
ear spectral properties of the optically driven evolution 
of the carrier subsystem which may decrease the overlap 
between these nonlinear spectra and the spectral features 
of the phonon reservoir. 

Instead of expressing the admissible class of pulses in 
terms of mathematically defined constraints, the pro- 
posed approach relies on physically transparent semi- 
quantitative conditions. By using pulses of finite dura- 
tion as building blocks one controls the frequency cut-off 
so that the spectra of the resulting pulses may be flex- 
ibly restricted to the frequency range of which one has 
reliable knowledge. In this way we guarantee that the 
result of the optimization is defined by the physical con- 
tent of the carrier-phonon model and not by its unphys- 
ical truncations. On the other hand, by optimizing the 
pulse amplitudes the spectral properties of the resulting 
evolution are chosen in such a way that resonances with 
phonon modes are minimized. 

Given a finite control time window, it is not possible 
to completely avoid the excitation of phonon modes and 
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Appendix 

UNIVERSALLY OPTIMAL 2-PULSE SEQUENCE 

In this appendix we show that the optimal amplitudes 
of two ultrashort pulses (in the limit of infinitely short 
pulses) depend only on the total rotation angle but not 
on the properties of the lattice reservoir. 

Let us denote the total rotation angle by a and the 
rotation angles performed by the first and second pulses 
(proportional to the corresponding pulse areas) by ai and 
ck2 = a — ai, respectively. Assuming the pulses to arrive 
a.t t — ±ttot/2 and inserting Eq. Ijlll) into Eqs. (|10|l and 
^ one finds 

• 2 ,0(1 witot \ / . 2 • ■ \ 

sm a + 2 I 1 — cos — - — I (sm ai — sm a sm ai ) . 
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Since 1 — cosut > 0, the minimum (with respect to ai) 
of the overlap with the nonnegative function R(uj)/uj'^ 
[Eq. lO] always corresponds to the minimum of the right- 
most bracket in the above formula. Therefore, irrespec- 
tive of the specific form of the optimal choice is 

1 . 

smai — — sma. 



For a = 7r/2 this corresponds to ai = tt/6 and a2 = 7r/3, 
so that the pulse area ratio is 1:2. This result holds, 
however, only in the leading order in the phonon coupling 
and is no longer valid in the general, non-perturbative 
case. 
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